A Landauer transport formulation with inelastic scattering 
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The Landauer transport formulation is generalized to the case of a dynamic scatterer with an 
arbitrary energy level structure, weakly coupled to a long ideal noninteracting wire. The two- 
terminal linear conductance of the device is expressed in terms of the (in general inelastic) scattering 
cross sections of the electrons off the scatterer, using a Fermi-liquid picture for the scattering states 
of the whole system. Assuming unitarity, i.e. the optical theorem, we obtain the same results from 
a Kubo linear response treatment and from a Keldysh technique approach. 
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Introduction. The scattering model due to Landauer 
^] has been extremely useful for studies of the electronic 
transport in mesoscopic systems, since it gives a clear 
physical picture of the processes involved (for reviews, 
see Refs.QandH. This formulation expresses the linear 
conductance ^ of a given finite sample (the scatterer) in 
terms of transmission probabilities for electrons injected 
to and collected from it, via ideally conducting leads. 
For noninteracting electrons, with a total transmission 
amplitude t and magnitude T — \t\'^, the two-terminal 
single-channel Landauer formula, including spin, reads 



2e^ 



T. 
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In this paper we generalize his result to include inelastic 
scattering by the sample. 

A general formal expression for the conductance of a 
scatterer imbedded in a long, ideal, "wire" was given in 
Ref. 4. It was proven to reduce to the Landauer for- 
mula at the zero temperature (T — > 0) limit. How- 
ever, at finite T it was not possible to produce such a 
proof, and the inelastic scattering channels, which open 
up then, were implicated. Later, the opinion that the 
Landauer formula is not applicable with interactions on 
the scatterer became prevalent (e.g. Ref. 0). That would 
mean that many of the usual mesoscopic effects would 
be severely modified by the inelastic scattering and in- 
teractions. Thus, the generalization (see also Ref. @) we 
propose is relevant. Indeed, the full scattering matrix 
of a local perturbation in an otherwise ideal conducting 
wire, should contain the necessary physical information 
for the conductance. The quasiparticle picture is nec- 
essary to approximately include the effects of electron- 
electron interactions on the leads mediated by virtual 
second-order processes on the scatterer. Within this pic- 
ture, multiparticle processes and non-Fermi-liquid behav- 
ior are avoided. 

The linear dc conductance of the sample, at finite 



temperatures, is simply expressed through its full single- 
quasiparticle scattering matrix (satisfying unitarity and 
time-reversal symmetries). The results are then con- 
firmed via the Kubo linear response approach , consid- 
ering the sample as a set of scattering states with equilib- 
rium populations. Finally, we use th e op tical theorem 
to express the linear Keldysh |j, Isl llll lul conductance 
in terms of the sample's transmission cross-sections, also 
confirming the generalized Landauer result. 
Scattering formulation for a one-dimensional 
(ID) conductor. We demonstrate our statement for 
the simple case of a scatterer (a "quantum dot"), taken 
as a finite "black box" with arbitrary interactions inside 
it, coupled weakly to ideal ID leads. We start by writing 
the scattering states, including inelastic processes. The 
states coming from the left {£) and from the right (r) 
belong to two independent Hilbert spaces, with no mu- 
tual Fermi blocking factors. For linear transport, the two 
chemical potentials of the {£) and (r) states, and /i,. 
respectively, are slightly different (/i^ < fJ-i — I^r = 
eV << fJ.e,r)- At equilibrium, both spaces are populated 
with the same temperature T and iii = fir — li- Without 
the coupling, the states \i) of the dot have discrete en- 
ergy levels, Ci, including the arbitrary interactions there. 
The leads have continuous spectra, characterized by a 
wavenumber fc, longitudinal velocity v = Tik/m and a ki- 
netic energy of the "incoming" electron E = h?k^ /2m. 
The states coming from the left lead (labelled by t) are 



eki) 




J), 
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far to the left of the dot and far on its right, Vj = 
hkj/m (> 0) is defined by ?i^fc|/2m = E — ej +ei, tji and 
rji are the transmission and reflection amplitudes from 
the left, exciting the dot from state \i) to state The 
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states coming from the right lead are 



j), and 
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far to the left of the scatterer and far on its right, with the 
transmission and reflection amplitudes t'j^ and r^j. Using 



T 

J- 1] 



the 



total right-going current in the states with an energy 
E va. a window dE, per spin direction, is 



ev TjiUidE 



(4) 



where rii = {2nfiv)^^ is the density of {iki} states per 
spin channel (which cancels the velocity factor, as usual). 

We next discuss the thermal occupations. The proba- 
bility to find the dot in state \i) is Pi, with X)i = 1- 
Consider a typical low-lying excited state of the whole 
system with the dot in state Within a Fermi liquid 
description the number of quasiparticles in this state is 
macroscopic. Thus the contribution of the dot energy 
to the energy of a single quasiparticle state is negligible. 
Hence their populations are given by the Fermi function 
of the electronic energy only. Thus, the populations of 
the scattering states | Iki) and | rki) are 



Peki,rki — Pif{E — fli,r), 



(5) 



where f{x) = l/{e^^ + 1) is the Fermi function and (3 = 
{kBT)~^. The total populations of incoming electrons 
with energy E are given by fi^r{E) = f{E - ^i^r)- 

The averaged right (left)-going current carried by the 
(, (r) states is given by 

li^r = Pjkihki — j dE'y^ PjkiTji; 

i ij 

= ^h^T.Prk.Tr. (6) 

ij 

For a finite, time-reversal symmetric dot, with real eigen- 
states, Tji — 7^'^, meaning that the cross-section for scat- 
tering of an electron from the dot, flipping the dot from 
state i to state j, has the same value irrespective of the 
direction from which the electron is coming. The total 
average net current is given by / = le^r — Ir^i- Using 
Eq. (jSJ, we find, to first order in eV = lie — iir, that the 
linear conductance, Q = limy^o I /V per spin is 



e 
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This is the two-terminal single-channel, finite-T Lan- 
dauer linear conductance formula including arbitrary in- 
teractions on the scatterer and inelastic scattering. 



Rederivation from the Kubo formula. Here, we con- 
firm our main result, Eq. (jTI, using the w — > limit of 
the Kubo linear response formulation 7], 



(8) 



The velocity power spectrum, Sv{uj) is given for Fermi 
liquid electrons in a unit volume and including spin by 

Sv{uj) EE ^ \Vmn?5{hL0 ~ Em + -E„) X 
m,n 

p^fiEm-^l)[^-fiEn~^l)], (9) 

where the matrix elements of the velocity operator v (in 
the wire direction) are taken between an initial state 
|m) — \£ki), or \rki) and a final one, \n) = \£k'i'), or 
\rk'i'). k is related to E ^ Em and k' to E' = E,,. The 
summation over m means summation over i, taking the 
sum of the i and the r terms and integrating over the 
energy E (which, as explained above, does not contain 
the dot energy), with the appropriate density of states, 
and similarly for the summation over n. Note that here 
the Fermi blocking factors, [1 — /(£-« ~ d,o appear, 
since real transitions are considered. The conductance is 
proportional to the net absorption rate from a classical 
field at frequency oj (12]. 

The transitions are of four types: (left-moving 
to left-moving), {i,r) (right-moving to left-moving), and 
similarly (r, r) and (r, £) . The sum of the two latter tran- 
sitions equals that for the two former ones, so we calcu- 
late only the former, and multiply the final result by two. 
From Eq. Q, 



{£ki\v\£k'i' 



ji' 1 
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where v' = hk'/m and we used the (column) unitar- 
ity properties of the full S-matrix, J^j tjitji' = ^i.i' " 
J^j^jiTji'. Likewise, {£ki\v\rk'i') = vWX;^- t*,r^,. Tak- 
ing the absolute values squared of these matrix elements, 
summing over i', (together with the density-of-states fac- 
tors required to convert the k- and fc'-sums into energy 
integrals) and using the (row) unitarity properties of the 



fuU S-matrix, X^i' = ^: 
Y,n.,ni>{\{£ki\v\£k'i')\'^ + \{£ki\v\rk'i')\ 



j,r , yields 



(11) 

Adding a similar contribution from the states which start 
in \rki), we put the resulting expression into Eq. (Q, and 
use f{E + hLu/2)[l - f{E - huj/2)] = f{E)[l - f{E)] + 
{hLu/2)f'{E) + 0{uj^), to obtain for smaU uj, 
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Eq. ^ now yields Eq. jT)) for Q{uj = 0), a,s required. 
Confirmation via the Keldysh technique. It has 

been claimed 0, Q that when the transmitted electron 
is subject to inelastic interactions on the scatterer, the 
transmitted current should be obtained using methods 
applicable to many-body problems. The most-frequently 
used one is the Keldysh technique. Here we calculate in- 
dependently the transmission of the system from the full 
scattering operator, and the current using the Keldysh 
technique, and then confirm that the two expressions so 
obtained are identical. Moreover, this equivalence is just 
a result of the optical theorem, namely, it follows from 
the assumed unitarity of the full scattering matrix. 

Considering a quantum dot with a single electronic 
level 

by the wavevector k and those on the right by p. Both 
leads are identical, except of being connected to reser- 
voirs of different chemical potentials, creates an elec- 
tron on the dot.) The current between the leads is given 
by the one flowing from the left, li — e{d/dt) J^ki'^lk'^^k) 
(where (...) is the quantum average), or by its counter- 
part —Ir- In the Keldysh technique these currents are 
expressed in terms of various Green functions, resulting 
in (see, e.g., Ref. 0) 

TfTr 



the coupling with the leads is J2k ^^kcl/^d + 
rpclpd+hc. (The band states on the left are denoted 



/ = - duj 



- fr 

where the partial widths. 



(cf -G;^)], (12) 
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are needed only in a narrow energy interval around the 
(average) Fermi level of the leads. Therefore, one may ne- 
glect the energy-dependence of these widths [3- G^'^{uj) 
are the thermal averages of the usual retarded and ad- 
vanced Green functions on the dot, which include all in- 
teractions. Despite the apparent similarity of Eq. H12|) to 
the Landauer formula, it has been claimed that there 
is no connection between the terms in the square brack- 
ets there and the transmission, once inelastic scattering 
is present. However, as we show below, these terms are 
in fact identical to the transmission. 

The transmission is derived from the full scattering 
operator of the system, S. Writing the total Hamilto- 
nian of the system as TC — Ho + Hi, where Hq is the 
Hamiltonian of the disconnected system (including the 
interactions on the dot) and Hi describes the coupling 



between the dot 
Hi), S = limt^oo 



uit) = i~- 



and the leads 
U{t), where 
t 

dt'e 



then (to all orders in 
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The transition probability per unit time between two 
states of the system, |m) and |n), is given by 



W„^„ = Ymi^({m\U^\n)j^{n\U\ 



m) + cc 



(15) 



When one specifies to the case in which both states |m) 
and \n) have no electron on the dot and differ only by 
one electron being on the left lead in the former and on 
the right lead in the latter, one finds |3. Iioj 



n JO 

At ^ „ _ \riR 



dTidT2dT^e n e R 



xGj(-r2 - n, -T2)G'^{0, -rg) + cc. 
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The Green functions appearing here depend on two time 
arguments (and not on their difference alone), since the 
thermodynamic average has not yet been performed. It 
is of crucial importance to realize that the Green func- 
tions appearing in the Keldysh technique result (|12|l for 
the (linear response) current are the Fourier transform of 
the thermodynamically averaged Green functions. On the 
other hand, the transition probability per unit time in- 
cludes a product of two Green functions. Obviously, the 
thermodynamic average of the product is not necessarily 
the product of the two averages. The relation between 
the two is supplied by the optical theorem. 

The transmission from a state of an electron with en- 
ergy E = huj on the left lead to a state of an elec- 
tron with energy E' = hto' on the right lead is 
T{lo ^ Lo') = p Wk^pSif - Lo)SC-i - uj')- Hence, the 
total transmission probability of an electron having ini- 
tial energy Hlu is T{u!) = J duj'T{uj — > oj'). For overall en- 
ergies near /i, T{uj) is identical to the quantity PiTji 
which appeared in Eq. jnj. With the definitions l|13|l . 
T(a;) becomes 



T{u;) 



dTidT2e^ 



'j{T2-ri) 



xGi{-Tl,0)G^{0,~T2), 



(17) 



and consequently, using the Landauer formulation (see 
also Ref. M) 



dcu 



(^Moj)^fr{cj))T{u;). (18) 



Next we show that the two expressions, Eqs. (I12f) and 
(|18|l . are identical, namely, that the terms in the square 
brackets in Eq. (^^J are given by T(cj), Eq. ^7^. To 
this end, we use Eq. (|14|l to write down explicitly the 
condition that S is unitary, '^^{n\S''\m) {m\S\l) — Sn.i- 
Specifying to the case in which the states \l) and \n) differ 
just by one electron taken from the left to the right lead, 
the unitarity condition becomes 
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Since Ti^r are independent of the energy, we have 

J2kp \yik,rp\'^e^'^* — '^Ti^r^t). Performing the thermo- 
dynamic average over Eq. (|19|l . taking into account that 
G^{h,t2)_{G^ih^) is defined for h > is {h < is), 
and that £'^'"^(^1,^2) depends only on {ti — <2), we find 



+ 2i- 



dte' 



:Giit-Tut)Gfit,t-T2] 



dTidT2e' 
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(20) 



Since the thermal average in the second term here is in- 
dependent of t, the i— integration yields 27r(5(efc — tp) and 
Eq. (1201) gives 



G^{ek)-Giiek)^-^n- 



TfTr 



(21) 



which shows that the square brackets in Eq. (|12|l are 
equal to T{uj), Eq. (|17|l . That is, the current using the 
Keldysh method, Eq. p2|) . and the same current by the 
Landauer-type formula, Eq. (fTH|l . are identical. We have 
confirmed this equivalence explicitly by considering the 
case in which the electrons on the dot interact with a 
phonon bath [l^ . 

Conclusions. Within the Landau framework (used here 
for singie-quasiparticle scattering states), we found that 
the Landauer-type picture for the linear transport is gen- 
eralizable to include inelastic scattering and interactions 
on the "dot" . We considered here a dot having, when it 
is closed, an arbitrary energy level structure, due to (pos- 
sibly strongly interacting) degrees of freedom (e.g. vibra- 
tions, spins) other than those of the conduction electrons. 
The results were confirmed with both the linear response 
and Keldysh approaches. Unitarity and time-reversal 
symmetries for the dot played crucial roles. However, 
the generalization for a dot imbedded in an Aharonov- 
Bohm interferometer is straightforward. So are those to 
multichannel and multiterminal situations. 
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